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Abstract 


This  paper  presents  a  simple  formula  for  calculating 
the  displacement  thickness,  momentum  thickness  and  physical 
thickness  of  a  laminar  boundary  layer  in  high-enthalpy  flows  in 
conic  nozzles  at  hypersonic  Mach  numbers.  The  inviscld  flow 
upstream  of  the  throat  is  assumed  to  be  in  chemical  and  vibra¬ 
tional  equilibrium,  and  downstream  of  the  throat  the  flow  is 
frozen.  The  present  expressions  for  6*,  0  and  6  are  given 
as  explicit  functions  of  the  geometric  dimensions  of  the  nozzle 
and  the  parameters  of  the  inviscid  flow  field  at  the  throat  and 
at  the  point  considered,  instead  of  an  integral  expression. 
Therefore,  it  is  particularly  suitable  for  the  design  of  this 
kind  of  nozzle.  Comparison  between  the  result  from  the  present 
formula  and  published  experimental  data  shows  that  the  agreement 
is  satisfactory. 


Symbols 

a  speed  of  sound 

A  cross-sectional  area  of  the  nozzle 
A  constant  in  equations  (10)  and  (17) 

B  constant  in  equations  (10)  and  (17) 

C  constant  in  equations  (17)  and  (18) 
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D  diameter  of  the  nozzle 

U  velocity  component  in  the  x-direction,  equation  (5) 
v  velocity  component  in  the  y-direction 

V  velocity  component  in  the  Y-direction,  equation  (5) 
x  coordinate  along  the  nozzle  wall 

y  coordinate  perpendicular  to  the  nozzle  wall 
X  transformed  coordinate,  equation  (2) 

Eq  constant,  equation  (14) 

/(m r,)  function,  equation  (15) 

asymptotic  function  of 
h  static  enthalpy 
H  total  enthalpy,  »“*  +  y 
Hfc  shape  factor,  h,,-**/*,, 

kQ  constant,  equation  (3) 
k,  ^  constant  equation  (14) 

1  shearing  stress  parameter, 

M  Mach  number 

n  pressure  gradient  correlation  coefficient, 
in  reference  [3] 

^  pressure  gradient  correlation  coefficient, 
equation  (7) 

N  momentum  parameter,  equation  (10) 

p  static  pressure 

P  Prandtl  number 
r 

r  radial  distance  of  a  point  inside  the  boundary 
layer  from  the  axis  nozzle 
r*  radius  of  the  nozzle  at  the  throat 
R  curvature  factor,  equation  (5) 
s  enthalpy  function,  equation  (8) 

T  static  temperature 
u  velocity  component  in  x-direction 

Y  transformed  coordinate 
z  compressibility  factor 


2 


*  semi-conical  angle  of  the  divergent  part  of  the  nozzle 
$  pressure  gradient  parameter,  reference  [31 

r  specific  heat  ratio 
r„  isentropic  exponent 

*  boundary  layer  thickness 

*•  displacement  thickness  of  the  boundary  layer 

*  integrated  value  in  Table  1 
I  momentum  thickness 

«  constant,  equation  (7) 

x  Sutherland's  viscosity  law  coefficient,  equation  (3) 
f>  viscosity  coefficient 

*  kinematic  viscosity  coefficient 
*•  constant,  equation  (7) 

<  area  ratio  of  the  nozzle  i 
0  density 

*»  semi-conical  angle  of  the  convergent  part  of  the  nozzle 

function, equation  (20) 
r,>  function,  equation  (21) 
r,)  function,  equation  (22) 
j.  stream  function 


Subscripts : 

o  stagnation  value  of  the  inviscid  flow 
e  parameter  at  the  outer  edge  of  the  boundary  layer 
eq  value  of  the  equilibrium  flow 
f  value  of  the  frozen  flow 

W  value  at  the  nozzle  wall 

tr(X,Y)  value  at  the  coordinates 
*  value  at  the  throat 


Whenever  the  coordinate  variables  x,  y,  X,  Y  appear  as 
subscripts,  a  partial  derivative  with  respect  to  that  variable  is 
implied . 
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I .  INTRODUCTION 

For  the  theoretical  analysis  of  a  boundary  layer  in  high- 
enthalpy  nozzle  flow  at  hypersonic  Mach  number,  only  approximation 
methods  are  applicable.  Until  now,  only  numerical  integration  for 
specified  stagnant  conditions,  geometrical  dimensions  and  wall 
temperature  could  be  done  by  means  of  sophisticated  computers  [1,2]. 
The  calculation  is  complicated,  and  the  results  obtained  by  dif¬ 
ferent  approximation  methods  are  usually  quite  different  from  each 
other  [3-8].  Because  of  the  dissimilarity  of  the  problem  and  the 
shortage  of  satisfactory  modifying  factors,  the  computed  results 
of  specific  nozzle  shape  and  stagnant  parameter  F  cannot  be  applied 
directly  to  another  nozzle  with  different  parameters.  On  the  other 
hand,  the  "low  density"  caused  by  the  "high  enthalpy"  and  "high 
Mach  number"  makes  the  boundary  layer  region  inside  the  nozzle 
even  larger  than  the  inviscid  flow  region.  This  means  that  a 
very  accurate  calculation  of  the  displacement  thickness  and  boundary 
layer  thickness  is  absolutely  required. 

It  was  not  until  1969  that  Petrie  [9]  did  some  experimental 
Investigations  on  high  enthalpy  laminar  boundary  layer  and  presented 
a  systematic  evaluation  to  decide  which  theoretical  method  Is  the 
most  desirable.  He  utilized  electro-beams  and  pressure  probes  to 
invest igage  the  boundary  layer  profile  inside  a  nozzle  flow  of  air 
with  stagnant  parameters  PQ  *  1  atmosphere,  Tq  *  5000  K,  and  dis¬ 
covered  that  all  theoretically  computed  values  for  6*  are  higher 
than  the  experimentally  monitored  value.  The  Cohen-Reshotko  [3] 
method  was  found  to  have  the  slightest  deviation,  sometimes  only 
45S  higher.  Since  then,  this  method  has  been  more  widely  applied 
for  the  computation  of  a  high  enthalpy  laminar  boundary  layer  with 
specified  nozzle  parameters  [14].  However,  it  is  very  difficult  to 
apply  the  method  in  the  practical  design  of  high  enthalpy  nozzles. 
This  is  because  the  integral  factor  In  the  equation  is  expressed 
in  terms  of  the  integral  form  of  the  inviscid  flow  parameters  of 
the  entire  nozzle  flow:  to  calculate  n,  the  outer  edge  parameter 
of  the  boundary  layer  must  be  integrated  numerically  from  the 
stagnation  point  at  the  nozzle  entrance  to  all  the  points  at  the 
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considered  cross-section  of  the  nozzle,  and,  in  the  course  of  the 
design,  all  calculations  have  to  be  repeated  along  with  more  than 
five  iterations  for  a  single  variation  of  any  of  the  nozzle  inlet 
parameters  pQ,  H  or  geometrical  parameters 
and  A/A*. 


Utilizing  the  fact  that  high-enthalpy  nozzle  flow  at 
hypersonic  Mach  numbers  is  generally  characterized  by  an  inviscid 
flow  in  equilibrium  upstream  of  the  throat  and  a  chemically  and 
vibrationally  frozen  flow  downstream  of  the  throat,  starting 
from  the  momentum  integral  equation  for  an  axial-symmetric 
boundary  layer,  modifying  the  definition  for  n  downstream  of 
the  throat,  applying  the  linearized  approximation  of  favorable 
pressure  gradient  universal  relation  similar  to  those  in  [3], 

[9],  [11]  and  [12],  employing  the  Mach  number  M  as  the  integra¬ 
tion  variable,  using  a  continuous  function,  partial  asymptotic 
approach  method  for  the  integral  function,  and  specifying  the 
condition  that  the  "cold  wall"  (being  a  high-enthalpy  ground 
simulation  wind  tunnel)  has  universality,  much  simpler  expressions 
for  n,  0,  <5*  and  <5  are  derived  in  this  paper. 

These  expressions  contain  explicit  functional  forms  of 
the  geometric  variables  of  the  nozzle,  parameters  at  the  throat 
and  at  the  cross  section.  Hence,  tedious  numerical  integration 
is  avoided. 


II.  GOVERNING  EQUATIONS  AND  BASIC  ASSUMPTIONS 


The  region  downstream  of  the  throat  is  axi-symmetric , 
steady  frozen  gas  flow.  The  boundary  layer  equations  are  written 
in  the  form  suggested  by  Probstein,  et  al  [13]: 


-—-(pur)  +  — —  (pvr)  —  0 
9*  9y 


du 

aT 


du 

"97 


pu 


dH 

d* 


dt>  Id/  d«  \ 

<•>*'*’*■  dy  vr  dy  ) 

dy  t  '  dy  '  Pr  dy  >  dy  L  \  Pr/  dy  i>  J 


(1) 
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in  which  H  *  h  +  5  . 
[15]  transformation: 


Introducing  the  Stewart son  [14]  and  Mangier 


(2) 


where  X  is  determined  by  the  viscosity  law: 


The  stream  function  is  defined  by: 


p  .  .  p 

-* - vr\  <p,  “  —  wr. 

Po  P 0 


The  following  relations  are  also  introduced: 


r  jly  -  ( 

.  vcosaV 

/  V 

Substituting  equations  (2)  to  (5)  into  (1)  and  simplifying, 
ux  +  vY  -  0 

0vx  +  Vyv  -  U.U.X  ^  ^  (RUy)r 

UHX+  VHy  -(1  -  J-)-25v0.e?^I^(2?UC;y)v  +  |?4«^^.(RHy)y 

\  Pr/  ?.•  p,  r„  Pr^;  A  r,  J 


Now  the  following  basic  assumptions  are  made: 

(a)  The  flow  inside  the  boundary  layer  is  also  chemically 
vibrationally  frozen.  According  to  experimental 
results  for  the  air  medium,  the  N£  vibrational  relaxa¬ 
tion  phenomenon  inside  the  boundary  layer  is  extremely 
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( 
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weak.  Hence,  it  can  be  assumed  that  the  chemical 
combination  is  even  weaker. 

b)  The  second  order  effect  of  the  lateral  curvature 
effect  is  neglected.  This  is  because  the  lateral 
curvature  effect  has  an  apparent  influence  on  the 
velocity  profile  and  the  friction  drag  but  very 
little  influence  on  6*.  Hence,  R  **  1  can  be  taken 
for  granted. 

c)  Pr  1.  Since  the  condition  of  a  strong,  favorable 
pressure  gradient  flow  at  a  hypersonic  Mach  number 
outside  an  isothermal  wall  is  involved  here,  Pr  =  1 
is  almost  valid.  Petrie  [9]  analyzed  the  velocity 
and  temperature  profiles  at  different  Pr  under  the 
condition  of  large  B  values.  The  observed  results 
by  means  of  an  electron  gun  show  that  the  assumption 
of  P  1  is  valid. 

d)  The  pressure  gradient  correlation  coefficient  n  is 
redefined  as: 


X2.  _ 

P*  Po  7" o 


Kva 


(7) 


e)  The  enthalpy  function  is  defined  as: 

,  _  *1f2  +  a1/(.v — 1)  _  j 
«;/2 +«;/(«•,- 1 ) 

and  it  is  assumed  that: 


t 


(8) 
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III.  THE  SIMPLIFICATION  AND  SOLUTION  FOR  THE  RELATION  FOR  n 

Under  the  above  basic  assumptions,  the  boundary  layer 
momentum  equation  can  be  integrated  along  the  Y-axis  from  0  to 
*  and  a  momentum  integral  relation  of  exactly  the  same  form  as 
those  given  by  the  references  [3],  [9],  [11]  and  [12]  is  ob¬ 
tained: 


The  universality  of  this  kind  of  differential  equation  for 
favorable  pressure  gradient  under  the  condition  of  incompres¬ 
sible  flow  was  first  proved  by  Thwaites  [11].  Later,  Kochln- 
Loytsyanskiy  [12]  and  Cohen-Reshotko  [3]  also  adopted  the  uni¬ 
versal  relation  of  this  form.  After  linear  approximation, 
equation  (9)  becomes: 


N(*w,  »)  ■»  +  2)  +  /]  •»  A  +  5n 


integrating,  the  expression  for  n  downstream  of  the  throat  is: 


-  -AU.xw;*  [[*  L 1‘r'dX  -  —  -  I 

IJ**  AUTxUT’1  i 


For  the  region  upstream  of  the  throat,  if  an  isentropic 
flow  of  a  perfect  gas  of  r«<l-400  is  assumed  to  be  equivalent 
to  the  isentropic  equilibrium  flow,  the  following  expression  is 
obtained: 


**  - -au:-*v:x  [[xV,-rfx-  — -  ;»  ] 


Modifying  (12)  by  the  factor  *  —  #  substituting  into 

Pi  Pi i  To 


TV.  : 


(11),  and  transforming  back  to  the  (x,y)  coordinates  plane, 
the  following  is  obtained: 


For  various  Yes  upstream  of  the  throat,  the  calcu¬ 
lated  values  of  the  first  integral  in  equation  (13)  are  shown  in 
Table  1.  When  ambient  pressure  pQ  is  varied  in  the  range  1-5 
atmospheres  and  H  in  the  range  2000-5000  kcal/kg,  the  values  of 
Yes  at  the  points  upstream  of  the  throat  vary  from  1.22-1.35 
correspondingly.  Based  on  Table  1,  the  variation  of  the  integral 
is  less  than  2.6%.  Also  froih  an  order  analysis,  the  second  in¬ 
tegral  in  equation  (13)  is  one  order  of  magnitude  higher  than  the 
first  integral  at  hypersonic  Mach  number.  Hence,  the  error  of  n 
due  to  variation  of  the  value  of  the  first  integral,  induced  by 
taking  different  Yoa  ,  is  generally  less  than  0.35$.  There- 

fore,  any  flow  of  perfect  gas  with  y  between  1.22  and  1.35 

©  s 

can  be  used  to  simulate  the  equilibrium  flow  and  calculate  the 
value  for  n.  Downstream  of  the  throat,  the  integration  variable 
x  of  the  second  Integral  in  equation  (14)  can  be  transformed  into 

>  by  replacing  the  isentropic  relation  for  a 
perfect  gas  with  the  isentropic  relation  for  the  chemically  and 
vibrationally  frozen  flow,  and  by  replacing  the  relation  for  the 
nozzle  flow  of  perfect  gas  with  the  relation  between  the  correspon¬ 
ding  area  ratio  A/A*  and  M  .  Transforming  E,  into  M  .  the  fol- 
lowing  is  obtained: 


(15) 


Taking  §  =  2.50,  the  first  term  in  equation  (15)  can  be  integrated 
for  any  .  For  the  second  integral  function: 


n  «-  1.400  is  taken  as  the  left  end  point,  and  different  am,s  Ar, 
are  expanded  and  asymptotically  approached  between  r>  ■■  1.400 
and  1.5.  In  doing  so,  the  asymptotic  function  for  /(M,,  r,)  at 
different  Mach  numbers  can  be  obtained: 


~~  -  0.167[  1  +  4.49(y,-  1.400)1  (1  <M.<2.5) 

“  °-°l6tl  +  6>25(r'  “  1-100)1  (2-5  <  M'  <  ,0) 
o  (M.  >  10> 


It  can  be  proved  that  when  1  <  M,  <2.5  ,  the  asymptotic 
error  (7 —  /)//<  2%.  For  2.5  <  M  <  5  ,  (7  —  /)//**  4.4 %  . 

fowever,  the  value  of  f  has  little  contribution  to  the 
integral  in  equation  (15  )  when  2.5  £  M  _<  5.  fence, 
it  can  be  shown  that  the  error  of  integration  is  still 
within  2%. 


Substituting  expression  (16)  into  equation  (15),  and  sub¬ 
stituting  the  n  value  corresponding  to  yes  =  1.25  along  with 
equation  (15)  into  (13),  and  applying  the  speed  of  sound  relation, 
the  differential  form  of  Euler's  equation  of  motion  and  the  energy 
equation  under  chemically  and  vibrationally  frozen  flow,  the  fol¬ 
lowing  can  be  obtained  after  simplification: 


n 


-A 


1  +  — •  M iffc' 


{J 

Mi  -  l  \ , 

<Pi(M,,  r()  —  »»  T/)| 


«in  a 
*in  <p 


(17) 
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where 


Xr„+» 

C"-(-rVr)'""'  (Sr.-l.25W, 


y„  -  1.25  B*t  C„  -  2.21) 


<PiUo)  —  0.0163  4-  0.4512  In  +  0.08|#-J 


ri) -  +  «y>* - (,  +  nfi  ^"P-| 

r,)  -  1.080  In  -----  --  4-  (0.883  -  0.75r, )f 2.5  - 

(‘W«/)  \  «;  / 

4 - - [2.5 _  M-J.Mr/-i.it*)1 

r,  — 1.166  J 

+  (0.124-  0.1 


Further  manipulation  gives  the  solution  for  momentum  thickness 
,  displacement  thickness  s#  and  boundary  layer  thickness 


e  -  -£*•  [ -^‘t-(i  +  — -  mjYT 

Tw  Uujdx  \  2  /  J 

5*  -  (9  4-  MKl  +  H.r)] 

d-d  \sje.r  4-  H-=-l-MKl  +  H„)] 


In  the  above  expressions,  and  8,r;e„  can  be  obtained  by 

the  similarity  profile  solution  of  the  Falkner-Skan  type  flow 
given  in  reference  [4];  «?,  is  determined  by  5*  ”*  j0  (*  —  jj  +  dY 

6„  is  determined  by  the  expression  8„  —  j"  ^  ft  —  dY  ; 

6<r  is  the  ij  — y  /5L±J.i£*.  value  corresponding  to  the  in- 

V  2  v#K 

compressible  plane  when  U/U  =  0.995,  and  n  is  the  value  of  n  at 
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which  NOw,m)  —  o  ,  corresponding  to  the  intersection  point 

of  8  *  1  and  s.  *  constant, 
w 

In  general,  water-cooled  copper  walls  are  used  for  high- 
enthalpy  flow  nozzles  in  simulation  wind  tunnels  on  the  ground. 

The  wall  temperature  satisfies  7V<400K.  •  When  the  stagnation 
temperature  T  is  between  2000  K  and  8000  K,  —  0.95  <  —  i  < 

O  l  o 

-0.80,  that  is,  the  following  relation  is  also  satisfied: 

- 1  <  ^  0.80  (26) 

Now  from  the  universal  relation  curve,  the  corresponding  value  of 
B  is  found  to  satisfy: 

2.36  <  fl  <  2.53 

From  equation  (11),  it  can  be  observed  that  the  B  index  inside  the 
integral  sign  Is  cancelled  by  the  B  index  preceding  the  Integral 
sign.  Hence,  It  is  regarded  that  the  variation  of  B  will  not 
cause  a  large  scale  variation  in  n.  Table  2  shows  the  difference 
in  n,  9  and  5*  caused  by  taking  B  as  2.5  and  2.3,  with  the 
same  nozzle  parameters.  Under  the  practical  condition  of  a  high- 
enthalpy  cold  wall,  as  discussed  in  this  paper,  the  value  of  B  Is 
closer  to  the  higher  extreme  of  the  Inequality.  Therefore,  the 
true  error  of  n,  9  and  6*  caused  by  taking  B  =  2.5  is  far 
smaller  than  those  in  Table  2.  Besides,  analyzing  the  data  from 
the  numerical  integration  of  the  similarity  solution  for  the 
momentum  integral  equation,  it  is  found  that  the  increase  in  B  also 
causes  an  Increase  of  A,  which  in  turn  reduces  the  variation  of  the 
value  of  n.  Hence,  it  is  reasonable  to  take  B  =  2.50. 

IV .  CONCLUSION 

In  the  present  effort,  the  inviscid  region  of  the  high- 
enthalpy  nozzle  flow  at  hypersonic  Mach  number  upstream  of  the 


throat  is  assumed  to  be  equilibrium  flow,  while  the  region  down¬ 
stream  of  the  throat  is  assumed  to  be  chemically  and  vibrationally 
frozen  flow.  The  flow  within  the  boundary  layer  is  also  treated 
as  chemically  and  vibrationally  frozen  flow,  on  the  basis  of  a 
monitored  data  analysis  of  the  N2  vibrational  temperature  relaxa¬ 
tion.  The  isentropic  relation  under  the  frozen  flow  condition  is 
utilized  for  the  outer  edge  of  the  boundary  layer.  The  modifi- 
cation  factor  is  introduced  to  modify  the  pres- 

sure  gradient  correlation  coefficient.  A  linear  approximation 
method  for  the  favorable  pressure  gradient  universal  relation 
is  applied,  and  the  cold  wall  condition  5  is  treated  as  single¬ 
valued.  Under  all  these  conditions  and  employing  order  of  mag¬ 
nitude  analysis,  integral  transformation,  dimensionless  Mach 
number  and  partial  asymptotic  approach  method  for  the  integral 
function,  the  explicit  functional  expressions  for  n,  6*  and 
6  are  derived. 

The  present  investigation  shows  that  when  pQ  is  between 
1-5  atmospheres  and  Hg  between  2000  and  5000  kcal/kg  or  varies 
within  an  even  greater  range,  the  calculation  of  the  pressure 
gradient  correlation  number  n  upstream  of  the  throat  can  be  done 
using  y  -  1.25,  or  the  "equivalence"  of  the  perfect  gas  flow 

^  5 

of  the  corresponding  isentropic  exponent.  Although  the  calcula¬ 
tion  of  n  downstream  of  the  throat  is  off  by  approximately  2% 
due  to  the  asymptotic  method  for  the  integral  function,  the 
results  at  the  nozzle  exit  calculated  by  the  formulae  presented 
here  are  still  close  to  the  experimental  data  after  modification 
of  the  isentropic  relation.  For  example,  for  air  at  TQ  =  5000  K, 

Po  =  1.16  atmosphere,  and  at  the  section  D/D*  =  9  downstream  of 
the  throat,  the  deviated  values  of  9  and  6*  are  corrected 
by  4. 2 %  with  a  modification  of  the  isentropic  relation  only. 

Table  3  lists  the  actual  samples  of  the  nozzle  calculation,  and 
Figure  2  shows  a  comparison  of  the  results  from  the  present  analysis 
the  results  from  other  theories  and  foreign  published  experimental 


Observing  the  results  in  Table  3,  it  is  apparent  that 
the  deviation  of  the  calculated  results  for  the  points  downstream 
of  the  throat  from  the  experimental  results  is  greater.  This  is 
quite  natural,  because  the  boundary  condition  has  switched  from  a 
non-slip  wall  condition  to  a  free  jet  flow  boundary  condition. 
This  portion  of  the  flow  field  should  be  calculated  by  a  non¬ 
similarity  method  of  free-shear  layer  [16].  The  methods  and  the 
results  given  here  are  no  longer  applicable. 


The  formulae  given  here  can  be  used  for  the  boundary  layer 
calculation  of  high-enthalpy  hypersonic  Mach  number,  cold  constant 
temperature  wall  conical  nozzles.  Since  common  integral  expressions 
are  avoided,  the  calculation  is  much  simplified.  Hence,  it  is 
especially  suitable  for  the  engineering  design  of  this  kind  of 
nozzle.  Common  high-enthalpy,  low-density  electric  arc  wind 
tunnel  nozzles  and  low-density,  low-temperature  plasma  flow 
nozzles  with  a  water-cooled  copper  wall  also  belong  to  this 
category . 


The  author  would  like  to  express  his  gratitude  toward  Mr. 
Yang-gwei  Bian  for  his  advice  and  assistance. 
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Figure  2.  The  growth  of  boundary  layer  displacement 
thickness  5*  :  q  Experimental  data,  1)  Pohlhausen 
energy  Integral;  2)  Pohlhausen  momentum  integral; 

3)  Beckwith-Cohen  method;  4)  Cohen-Reshotko  method; 
5)  modified  method,  reference  [1];  and  6)  equations 
17  -  25  given  here. 
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1)  Table  3  Comparison  of  Final  Results  With  an  Exclusion  thickness 

of  $• 

2)  Exclusion  thickness  (cm);  3)  Experimentally  determined  values 

4)  Results  with  an  IBM  7094  computer  using  the  Cohen-Resnorko  method 

5)  Formula  (16),  and  so  on,  as  found  in  this  article; 
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Abstract 

This  paper  presents  a  simple  formula  for  calculating  the  displacement  thickness, 
momentum  thickness  and  physical  thickness  of  laminar  boundary  layer  in  high  enthalpy 
flows  in  conic  nozzles  at  hypersonic  Mach  numbers.  The  inviscid  flow  upstream  of 
the  throat  is  assumed  to  be  in  chemical  and  vibrational  equilibrium  and  downstream 
of  the  throat  the  flow  is  frozen.  The  present  expressions  for  d*,  0  and  d  are  given 
as  explicit  functions  of  the  geometric  dimensions  of  the  nozzle  and  the  parameters  of 
the  invicid  flow  field  at  the  throat  and  at  the  point  considered,  instead  of  an  integral 
expression.  Therefore  it  is  particularly  suitable  for  the  design  of  this  kind  of  nozzle. 
Comparison  between  the  result  from  the  present  formula  and  published  experimental 
data  shows  that  the  agreement  is  satisfactory. 


